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Abstract: Detached eddy simulation (DES) based on the v 2 -f RANS model is pro- 
posed. This RANS model incorporates the anisotropy of near- wall turbulence which is 
absent in other RANS models commonly used in the DES community. In LES mode, 
the proposed DES formulation reduces to a transport equation for the subgrid-scale 
kinetic energy. The constant, CdeSi required by this model was calibrated by simulat- 
ing isotropic turbulence. In the final paper, DES simulations of canonical separated 
flows will be presented. 
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1 Introduction 

Detached-eddy simulation (DES) is a hybrid RANS /LES approach that performs RANS in attached regions and LES in detached 
regions using a single model. In the LES region, the length scale of the RANS model is set proportional to the grid size A. The 
RANS model thereby becomes an LES model. 

Although the Spalart-Allmaras (SA) model has been widely used for DES [1], its near- wall damping does not distinguish 
between velocity components. In contrast, the v 2 -f formulation [2] models the suppression of wall normal velocity fluctuation 
caused by non-local pressure-strain effects. This anisotropy has been shown to improve prediction of separation and reattachment 
[3, 4]. In addition, the SA model uses the minimum distance to the wall as the turbulence length scale which is not necessarily 
accurate at or near separation. The v 2 -f model, on the other hand, computes a length scale based on flow properties (i.e. the 
kinetic energy k, the dissipation rate e, and the kinematic viscosity v). Finally, in the LES region, the present treatment gives 
a transport equation for subgrid-scale (sgs) kinetic energy which is less empirical than the sgs viscosity transport equation used 
in SA based DES. 

2 The Proposed v 2 -f DES Model and Preliminary Results 

The v 2 -f formulation [2] has transport equations for k, e and v 2 and an elliptic relaxation equation for a function / analogous 
to the Reynolds stress redistribution term. These equations contain length and time scales L and T. The choice between RANS 
and LES modes is made by setting these scales appropriately. When the grid is fine enough to capture large turbulent eddies, 
i.e. when CdesA < fc 3 / 2 /e, the LES mode is selected. Otherwise, the RANS mode is selected. As suggested by Spalart et al. [5], 
the DES coefficient Cdes is chosen to match the correct energy spectrum in isotropic turbulence. In LES mode, the time scale is 
obtained as T = CdesA/a/A: where k now represents the sgs kinetic energy. Although the k equation is sufficient for purely LES 
simulations, the other three equations are required for the RANS mode. The v 2 -f model [2] is modified so that the entire set 
of the equations reduces to a transport equation for the sgs k in LES mode. In particular, coefficients in the elliptic relaxation 
equation are modified so that v 2 becomes (2/3) k in LES for isotropic turbulence. The final paper will include the formulation 
in detail. 

To validate the modified coefficients, fully developed turbulent channel flows with friction Reynolds numbers Re T up to 2000 
are simulated in RANS mode. Figure 1 shows that this modification maintains the performance of the unmodified RANS model 
with better prediction of v 2 near the center. Next, to check the LES mode and calibrate CdeSi decaying isotropic turbulence 
is considered. The initial field uses 512 3 DNS data [6] at Re\ = 105 filtered down to 32 s . Turbulence quantities in the DES 
model are initialized by running for some time with a frozen initial velocity field. The flow is allowed to run about eleven initial 
large-eddy turnover times at which point Re\ = 61. The v 2 -f DES simulation with Cdes = 0.8 agrees very well with DNS [6] 
and experiments [7] (figure 2(a)) up to the cut-off wave number. Because it allows a higher Reynold number and a wider inertial 
range, forced isotropic turbulence at Re\ = 98 is also considered using the stochastic forcing of Eswaran and Pope [8]. The v 2 -f 
DES gives almost exactly the energy spectrum produced by the dynamic Smagorinsky model (figure 2(b)). 

3 Conclusions and Future Work 

DES based on the v 2 -f model was proposed and validated for two cases: purely RANS and isotropic LES. These preliminary 
results are encouraging. The final paper will include DES simulations of the backward-racing step flow [9], and v 2 -f based DES 
will be compared with SA based DES. 
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Figure 1: Profiles of U + and v +2 of DNS(dots) and RANS with the unmodified [2] ( — ) and the modified v2f model in 
turbulent channel flows. The plus superscript indicates normalization with wall units. 




(a) Decaying isotropic turbulence ( Re\ = 61). 


(b) Forced isotropic turbulence ( Re\ = 98). 


Figure 2: Energy spectrum of isotropic turbulence. The proposed DES model (V2F-DES) is simulated with Cdes = 0.8. 
Kolmogorov -5/3 spectrum indicates E(k) = Ce 2//3 « -5 / 3 with the Kolmogorov constant C = 1.5. DynSmag indicates LES with 
the dynamic Smagorinsky model. 


References 

[1] Spalart, P. R., “Detached-Eddy Simulation,” Annual Review of Fluid Mechanics , Vol. 41, 2009, pp. 181-202. 

[2] Lien, F. S., Kalitzin, G., and Durbin, P. A., “RANS modeling for compressible and transitional flows,” Proceedings of the 
Summer Program, Center for Turbulence Research, Stanford, 1998. 

[3] Durbin, P. A. and Pettersson Reif, B. A., Statistical Theory and Modeling for Turbulent Flows, chap. 7.3 and 7.4, Wiley, 
2001 . 

[4] Iaccarino, G., Marongiu, C., Catalano, P., and Amato, M., “RANS simulation of the separated flow over a bump with active 
control,” Tech, rep., Center for Turbulence Research, Stanford University, 2003. 

[5] Spalart, P. R., Jou, W. H., Strelets, M., and Allmaras, S. R., Comments on the feasibility of LES for wings and on a hybrid 
RANS/LES approach, Greyden Press, 1997. 

[6] Wray, A. A., “HOM02: Decaying Isotropic Turbulence,” A Selection of Test Cases for the Validation of Large-Eddy Simula- 
tions of Turbulent Flows, AGARD Advisory Report No 345, 1997. 

[7] Comte-Bellot, G. and Corrsin, S., “Simple Eulerian time correlation of full- and narrow-band velocity signals in grid-generated 
‘Isotropic’ turbulence,” J. Fluid Mech., Vol. 48, 1971, pp. 273-337. 

[8] Eswaran, V. and Pope, S. B., “An examination of forcing in direct numerical simulations of turbulence,” Computers & Fluids, 
Vol. 16, No. 3, 1988, pp. 257-278. 

[9] Le, H., Moin, P., and Kim, J., “Direct numerical simulation of turbulent flow over a backward-facing step,” J. Fluid Mech., 
Vol. 330, 1997, pp. 349-374. 


